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It is pointed out that the unmagnetized inhomogeneous plasmas can support a low 
frequency electromagnetic ion wave as a normal mode like Alfven wave of magnetized 
plasmas. But this is a coupled mode produced by the mixing of longitudinal and 
transverse components of perturbed electric field due to density inhomogeneity. The 
ion acoustic wave does not remain electrostatic in non-uniform plasmas. On the 
other hand, a low frequency electrostatic wave can also exist in the pure electron 
plasmas. But the magnetic field fluctuations in both electron as well as in electron- ion 
plasmas are coupled with the electrostatic perturbations in unmagnetized case. The 
main instability condition for these low frequency electrostatic and electromagnetic 
modes is the same |fc n < kt (where K n and kt are inverse of the scale lengths of 
density and electron temperature, respectively). 

PACS numbers: 

I. INTRODUCTION 

A low frequency electromagnetic wave has been investigated as a fundamental normal 
mode of unmagnetized inhomogeneous plasmas. The ion acoustic wave (IAW) is a well- 
known low frequency mode of unmagnetized plasmas but it is purely electrostatic. Even for 
lightest ions of hydrogen the electron to ion mass ratio is very small — ~ 10~ 3 << 1. In 
the limit m e — > 0, the electrons are assumed to follow the Boltzmann density distribution in 
the electrostatic field of the IAW. 

In fact the electron inertia can play a vital role in producing a low frequency electromagnetic 
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wave. More than a decade ago [1], it was proposed that an electromagnetic wave having fre- 
quency near IAW is a normal mode of unmagnetized plasmas which can be responsible for 
magnetic field generation if it becomes unstable. In the derivation of its linear dispersion 
relation the electron inertia was not ignored while the displacement current was neglected. 
It was shown that the compressibility and vorticity can couple due to density inhomogeneity 
and hence a low frequency electromagnetic wave can be produced. Basically electrostatic 
IAW and magnetostatic mode [2] cooperate with each other to develop such a wave. The 
electron temperature perturbation was not taken into account and the steady state was as- 
sumed to be maintained by external mechanisms. The theory was applied to explain the 
magnetic field generation in laser plasmas. But the longitudinal and transverse characters of 
electric field decouple if the quasi-neutrality is used pp. Then both the high frequency and 
low frequency electromagnetic instabilities were also investigated in unmagnetized plasmas 

®- 

The interest in the investigation of low frequency magnetic fluctuations in unmagnetized 
plasmas was initiated after the first experimental observation of magnetic field generated in 
a laser induced plasma |4J. Then more experiments were performed |5j E] on these lines. 
Several theoretical models were presented to explain the magnetic field generation in laser 
produced plasmas [THT3]. 

The magnetic electron drift vortex (MEDV) mode was proposed as a pure transverse linear 
mode which can exist because of the electron temperature fluctuations in unmagnetized in- 
homogeneous electron plasmas [T3] . This mode can become unstable [2] if the equilibrium 
electron temperature gradient is parallel to the density gradient and it is maintained by ex- 
ternal effects. The basic MEDV mode is believed to exist in an electron plasma with smooth 
gradients. However, it has also been shown that the instability of MEDV mode can arise 
when the density profile is represented by a single step connecting two regions of nonzero 
density and temperature gradients [13]. The instabilities of magnetic and acoustic waves 
driven by perturbed baroclinic vector have also been investigated in a pure electron plasma 
|16j . Recently the nonlinear evolution of two dimentsional MEDV modes has been studied 
using computer simulation [TTj . The spontaneous magnetic field generation and formation 
of nonlinear structures has been discussed in detail. Most of the mechanisms proposed to 
explain magnetic fluctuations in initially unmagnetized plasmas are based on electron mag- 
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netohydrodynamics (EMHD) which has been discussed in detail in Refs. [T5J H5] . Some 
weaknesses and contradictions of EMHD model were pointed out several years ago [20] ■ The 
long-lived and slowly propagating nonlinear whistler structures (NLWS) or whistler sphero- 
maks (WSPS) have also been studied [21] using EMHD equations. Such structures have 
been observed in magnetized laboratory plasmas [221 123] . 

In MEDV mode the divergence of electric field is assumed to be zero (V.Ei) = while the di- 
vergence of electron velocity is non-zero (V.v el ^ 0). Furthermore the ions are treated to be 
stationary. The frequency u> of the MEDV mode is assumed to lie in between the ion plasma 
frequency and electron plasma frequency i.e. u pi « uj « u pe where u P j = J, for 

j=e,i and c is speed of light while k is the wave vector. These restrictions and assumptions 
are indeed very strict [20] and are not fulfilled in general. 

Several authors have considered the role of ion dynamics in the magnetic instabilities in un- 
magnetized plasmas. The effects of ion dynamics on MEDV mode have been investigated in 
the frame work of local approximation [21] . The coupling of high frequency electromagnetic 
wave with low frequency ion acoustic wave has been discussed in certain limits [3] ■ The cou- 
pling of magnetic fluctuations with ion acoustic wave has been discussed in a plasma with the 
steady state given as Vp e o = [25] without considering electron temperature perturbation. 
But the group velocity turns out to be negative in this treatment. 

It is important to find out some electromagnetic mode taking into account the ion dynamics 
using minimum approximations so that the strict restrictions on the frequency and wave- 
length of the perturbation are relaxed. It is better if the only required condition on frequency 
becomes u « u pe , ck. If we assume a steady state as Vp e o = 0, then the temperature gra- 
dient becomes anti-parallel to density. But laser and astrophysical plasmas are open systems 
and many external mechanisms can maintain a study state with parallel density and tem- 
perature gradients. For example, in stellar cores, both the density and temperature increase 
towards centre of the star due to fusion and star is held intact because of gravity. Therefore 
both the cases of parallel and anti-parallel gradients should be discussed. The electron ther- 
mal fluctuations can produce electromagnetic wave as was proposed many decades ago [T3] 
but the assumptions used in this work are very restrictive. 

It is not necessary to assume a pure transverse perturbation in an electron plasma. Rather 
the perturbation can be partially transverse and partially longitudinal and hence the electron 
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density perturbation may not be neglected. The frequency of such a wave in electron plasma 
turns out to be near (X e k y )v te K n where v te = f^M)^ is the electron thermal speed, A e = 
is the electron skin depth and the inverse of density gradient scale length 

L n — In the local approximation we need to have k„ << k. Furthermore the condition 

« Vf e k 2 can be satisfied if ^ << \ 2 De ky. Therefore, generally we may have v te ^ n < &pi 

i 

and hence one cannot neglect ion dynamics. Moreover, Vt e ^ n — c s k (where c s = 0^) 2 is i° n 
sound speed) is also possible. Therefore, it is expected that the electron thermal fluctuations 
can couple with IAW to produce stable and unstable low frequency electromagnetic waves 
in unmagnetized plasmas. It will be shown that the dispersion relation of ion acoustic wave 
is modified in the inhomogeneous plasma due to the coupling of electrostatic and magnetic 
fluctuations. A low frequency electrostatic mode can do also exist in a non-uniform pure 
electron plasma. This mode can couple with the IAW in electron-ion plasma. Several low 
frequency electrostatic and electromagnetic waves of inhomogeneous unmagnetized electron 
and electron-ion plasmas are investigated. Interestingly the main instability condition for 
these modes is the same. 



II. LOW FREQUENCY ELECTROMAGNETIC WAVES IN ELECTRON 

PLASMAS 

Let us consider the electron plasma in the background of stationary ions. First we discuss 
the dispersion relation of pure transverse MEDV mode [131 E]- Then we show that the 
compressibility cannot be neglected. Finally linear dispersion relations of low frequency 
electrostatic and electromagnetic perturbations are obtained in an electron plasma. The set 
of equations for MEDV mode in the linear limit can be written as, 

m e n d t v el = -en Ei - Vp e i (1) 
Air 

V x Bx = — Jj (2) 

c 

Ji = -en v el (3) 

VxE 1 = -i^Bi (4) 

Since pi = n T e i therefore energy equation becomes, 

3 3 

-n d t T e i + -n (v e i.V)T e o = -p V.v e i (5) 
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Curl of (1) gives, 

dtiV x v el ) = — d t Bi + — (Vn ) x VT el (6) 
m e c m e n 



Equations (2) and (3) yield, 



and hence 



Vel = -T^ V X B l ( 7 ) 



V x v e i = V 2 Bi (8) 

where Vrio x (V x Bi) = due to the assumption kJLVno-LBi. Equation (8) predicts 
E x = Eik. while Vn = x^p, V = (0, ik y , 0) and Bi = Biz have been chosen. 
Equations (6) and (8) yield, 

(1 + X 2 e k 2 )d t Bi = --(Vn x VTei) (9) 
en 

where A e = Equation (7) yields, 



V.Vel = — .(V X Bi) (10) 

and therefore one obtains, 

2 c 

T e i = k v K n Bi (11) 

Then (9) and (11) give the linear dispersion relation of MEDV mode as, 

uj 2 = \Co^fv 2 Te kl (12) 

where Co = j^lp an d v te — (T e /m e )^. The geometry of MEDV mode in cartesian co- 
ordinates is shown in Fig. 1. 
If VT e0 7^ is assumed, then (12) becomes, 



, ,2 _ r K n 

u — Co — 

Ky 



VtJT (13) 



where k t = l^-^^l and VT e0 = H-x^p 1 has been used. If (|/c n < n T ), then the mode 



I T e0 dx I ""^ v eu 1 ~ da; 

becomes unstable 



It can be noticed that V.Ei = and V.v e i ^ have been assumed in the above treatment 
and it does not seem to be very convincing. If we write (1) in x and y components, we can 
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FIG. 1: The MEDV mode geometry shows that it is a pure transverse mode while the electron 
equation of motion indicates E y \ ^ and hence k.Ei ^ 0. 

observe that the y-component of Ei should not be considered as zero due to d y p e \ ^ and 
hence the wave may be partially transverse and partially longitudinal. Equation (1) yields, 

e 2 / r el , . % 



and 



d t Vexl = E xl - v Te < K n — + k t — \ (14) 

m e I ±q no 



V V ^0 ^0 

It is obvious from (15), that V.v el ^ which implies E yi = —d y ^>\ 7^ 0. It is clear from Fig. 
1 that MEDV mode should be revisited including longitudinal effects. 
The curl of (2) yields a relation between E lx and E iy foru« u pe as, 

E u = -~(iE ly ) (16) 

(X Ky 
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where a = (1 + A 2 &; 2 ). 

Using equation of motion (1) instead of (10) in equation (5), we find, 



W. 



2±el 

T 



Wei 



e / , „ 3 „ 
n m e V 2 



(17) 



where W 2 = §u; 2 - A* (l - § and r 2 = (*r-y)*r . 

The continuity equation yields, 



-2 n el 

n 



V 2 k 2 



W 2 



3 k t v 2 e k 2 



2 /;„2- 

;</- 



-kyEly 



1 + 



(18) 



where L 2 = [u 2 - v 2 e k 2 (l - K 2 Jk 2 ) - ^(1 - (l - §f - T 2 ) }. The Poisson equa- 

(19) 



tion 



V.Ei = -47re 



n 



can be written as, 



ik y E ly [L 2 W 2 - u 2 pe {^ 2 + v 2 e k 2 y (l - K 2 Jk 2 y )}] 



(K n E lx )u 



pe 



(20) 



Equations (16) and (20) yield a linear dispersion relation in the limit u 2 « u 2 e as, 

a [L 2 W 2 - u 2 pe {W 2 + v 2 e k 2 y (1 - K 2 Jk 2 ) }] 



-U 2 )e {K n /kyf 



w 2 + 



3 kt 

\ 

2 K n 



v :2 



This equation can be simplified as, 



iW 2 = -~V 2 e K 2 n 



1 _ 



where 



{^ + \ 2 De k 2 y ) + 2 -\ 2 De k 2 y )a-K 2 Jk 2 y 



(21) 



(22) 



1 \ y g 

The second term in right hand side of (22) will disappear if E ly = is assumed. The first 
term and the factor \ 2 k 2 in second term are the contributions of transverse components 
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E\ x . Thus one can not obtain MEDV-mode dispersion relation for E\ y = from (22). It is 
interesting to note that equation (22) can be simplified to obtain, 

= ^lk 2 y (vlK n ) (l - ^Kr/K^j (23) 

which looks very similar to MEDV mode dispersion relation. 

However, the instability condition for this electromagnetic mode is the same as was for 
MEDV mode that is 

2 

-K n < k t (24) 

The low frequency mode (23) is partially transverse and partially longitudinal. 

If Ei x = is assumed, then equation (21) yields a low frequency electrostatic wave in a 

non-uniform unmagnetized plasma with the dispersion relation, 



2 _ V te r -n \ :j -"a ■"«■/ / 2 rx 
[(1 + Xbe k y) + l X De k f\ 



The instability condition remains the same (24). This indicates that the instability predicted 
by the so called pure transverse MEDV mode is always coupled with electrostatic perturba- 
tions in electron plasmas. 

In electron-ion plasmas, these modes of equations (23) and (25) can couple with the ion 
acoustic wave. 

III. IAW AND MAGNETOSTATIC MODE 

Here we shall show that transverse magnetostatic mode j2] which is obtained in the 
limit uP « Up e can couple with IAW in a nonuniform plasma [TJ. Therefore, both ions 
and electrons are considered to be dynamic while the electron temperature perturbation is 
neglected. The ions are assumed to be cold for simplicity and therefore equation of motion 
becomes, 

d t v a = —Ex (26) 
rrii 



The continuity equation yields, 



n rriiUj 2 



[k u E 1x + ik y E ly ) (27) 
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If quasi-neutrality is used then the transverse component E\ x and longitudinal component 
become uncoupled PQ. Therefore, we assume X 2 De k 2 7^ and use Poisson equation which 
in the limit u 2 « u pe becomes, 

[-v 2 te ky - u 2 pt (u 2 - v 2 e k 2 y ) - u 2 pe u 2 } ik y E ly ~ [u 2 pi (uj 2 - v 2 e k 2 y ) + u 2 pe u 2 ]K n E lx (28) 

Note that the term v 2 e k 2 is not ignored compared to u pe to couple E lx and Ei y . 
Equations (16) and (28) give a linear dispersion relation as, 

c^(a-^f) 
to = ^— 

where c 2 = ^% and b = (1 + \ 2 De k 2 )- We have to use Poisson equation to obtain a quadratic 
equation in u while (2) implies n el ~ tin. The equation (29) is the same as equation (21) 
of Ref. P] where two small terms in the denominator are missing. Note that X 2 De < \ 2 e and 
if quasi neutrality is used due to Ampere's law, then (29) yields the basic electrostatic IAW 
dispersion relation u 2 = c 2 k y . If the displacement current is retained and Poisson equation 
is used without using u 2 , u pi « u pe , then one obtains a dispersion relation of coupled three 
waves; ion acoustic wave, electron plasma wave and high frequency transverse wave [3]. 
Actually the contribution of displacement current in the curl of Maxwell's equation has 
been neglected for uj 2 « u pe ,c 2 k 2 . This should not mean that the electrostatic part of 
current is also divergence free, in our opinion. In the divergence part of Maxwell's equation 
lo 2 « u pe is used but v 2 e k 2 term is assumed to be important. It may be mentioned that 
in this treatment, Ampere's law does not imply quasi-neutrality necessarily. We need a 
coupling of divergence part and curl part of the current and for this we need to assume 
^ < X 2 nfi k 2 in the limit u 2 « u 2 , c 2 k 2 . 

rrii ue y pel 

IV. LOW FREQUENCY ELECTROMAGNETIC ION WAVES 

Now we present a simple but interesting theoretical model for low frequency electromag- 
netic waves assuming ions to be cold. The electrostatic waves will also be considered and 
it will be shown that magnetic field perturbation is coupled with the dominant electrostatic 
field. In the low frequency limit \dt\ << u pe ,ck, the electrons are commonly assumed to be 
inertial-less, i.e. — — > 0. Then the longitudinal and transverse components of electric field 



(29) 
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decouple. For IAW it is assumed that E = — Vy? while electron equation of motion yields 
Boltzmann density distribution as 

^ ~ e-f (30) 
n 

Then the fundamental low frequency mode of the plasma turns out to be the ion acoustic 
wave with linear dispersion relation 



u 2 s = c\k 2 y 



(31) 



in the quasi-neutrality limit. If dispersion effects are included, then instead of (31) one 
obtains, 

r 2 k 2 

(32) 



In the limit 1 << \ 2 De k 2 equation (32) gives ion plasma oscillations u 2 = u 2 ^ It is important 



<< 



to note that in the presence of inhomogeneity, a new scale ^ is added to the system. If : 

j , then longitudinal and transverse components of electric field can couple to generate 
low frequency electromagnetic waves. The divergence and curl of (1) give,respectively, 



d t V.(n v e i) = -— noV.Ei - — Vn .E! - — (V.Vp e i) 

m e m e m e 



and 



d t (y X V e i) + (K n x d t v el ) = K n X E x 



m P 



V x Ei 



(33) 



(34) 



where K n = |^^f| and Vn = +x|^£f| has been assumed. If initially electric field was purely 



dx 



electrostatic i.e. Ei = — V</?i, then it will develop a rotating part as well if Vn x Ei ^ 
as is indicated by the right hand side (RHS) of (34). 
The Poisson equation in this case is, 

V. El = 47 rn e - 

Using (18) and (27), the above equation can be written as, 



(35) 



ik v E 



L 2 W 2 u 2 - L 2 W 2 oo 2 pt - co 2 W 2 + v 2 te k 



— K nE\ x 

Then (16) and (36) yield, 



3K T 2 2 



i _ s. 

U2 



K y 



(36) 



a[L 2 W 2 -u 2 pe (W 2 + v 2 e k 2 y )(l- K 2 n /k 2 y )] 
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h.1 

K y 



pe 



Wn + 



3 k,t 



2 K n 



(1 - K°/*2) 



, ,2 

H ^— W j 



a — 



i-2 



In the limit cu 2 ,Wpj << io 2 e , (37) gives a linear dispersion relation, 



pe 



1 

Ho 



+ [a K>n/ky) C l K y 



K 



(37) 



(38) 



3 

If ion dynamics is ignored then (38) reduces to (23). The above equation shows a coupling of 
ion acoustic wave with the electromagnetic fluctuations in nonuniform unmagnetized plas- 
mas. 

In the electrostatic limit (E lx = 0), equation (38) becomes, 



1 

H~i 



2 2 
V te K n 



- K T /K n +C 2 s k 2 < - + 



(39) 



where H± — {(1 + A|, e A; 2 ) + §A|> e A^}. For stationary ions, (39) reduces to (25). 

It is important to note that the main instability conditions for electrostatic ion acoustic wave 

(39) and low frequency electromagnetic wave (38) is again the same (24). 

In electron-ion plasma, the electromagnetic wave dispersion relation in the quasi-neutrality 

limit can be written as, 



u> 



X^ky 



2 2 
; V te K n 



{a -Kim 

The wave geometry is shown in Fig. 2. 

The instability can occur if (24) is satisfied along with 
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(40) 



\ 2 k 2 

„2>,2 < 

2 s y (a- Kl/kl) 



2 2 



(41) 



V. DISCUSSION 

The theoretical model presented here shows that several electrostatic and electromagnetic 
low frequency waves can exist in un-magnetized electron as well as electron-ion plasmas. In- 
terestingly the main instability condition for these modes is the same |« n < kt where n n and 
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z 

t B, 




FIG. 2: The simplest possible geometry of electromagnetic ion wave is shown. This low frequency 
wave is partially longitudinal and partially transverse. 

k,t are the inverse of density and electron temperature gradient scale lengths, respectively. 
This indicates that the magnetic field fluctuations are always coupled with the electrostatic 
perturbations. 

It is well-known that the ion acoustic wave is a fundamental low frequency electrostatic mode 
of un-magnetized plasmas. Here we have found that the inhomogeneous electron plasma can 
also support a low frequency electrostatic mode. This electrostatic mode can couple with 
the magnetic field perturbations to give rise to a partially electrostatic and partially trans- 
verse wave. Therefore instead of the so called magnetic electron drift vortex (MEDV) mode 
there exists a low frequency electromagnetic wave having both the contributions of longi- 
tudinal and transverse electric field components as has been shown in equations (22). The 
dispersion relation of MEDV mode is very similar to the electromagnetic wave discussed in 



13 



this investigation as can be seen in equation (23). But the important point to note is that 
this dispersion relation appears after a cancellation of two terms. One of these terms is a 
part of longitudinal electric field E\ y . Moreover, if transverse electric field component E\ x 
is neglected one obtains a pure electrostatic wave of equation (25). The frequency range of 
both the modes is very close to each other. The instability conditions are almost the same. 
This fact strengths the view point that magnetic fluctuations are coupled with the dominant 
electrostatic fields. 

Any initial electrostatic field perturbation can produce its transverse component in the pres- 
ence of density gradient [TJ [3] . This phenomenon can cause a coupling of ion acoustic wave 
(IAW) with the low frequency transverse magnetostatic mode. This coupled mode has al- 
ready been investigated more than a decade ago pQ. But it can exist in a relatively shorter 
wavelength range for ^ < \ 2 De k 2 . In the quasi-neutrality limit, the IAW and magnetostatic 
modes decouple. 

In case of the electron-ion plasmas, the IAW does not remain electrostatic in inhomogeneous 
plasmas [IJESIEE]- The electromagnetic mode discussed for the case of pure electron plasma 
can couple with ion acoustic wave as shown in equation (40). Similar to the electrostatic 
IAW, this low frequency electromagnetic ion wave can exist even in the quasi-neutrality 
limit. The electrostatic and electromagnetic waves discussed here can become unstable if 
the density and temperature gradients are parallel to each other which can be the case in 
laser plasmas similar to stellar cores. 

These low frequency waves can be the intrinsic source of magnetic fields in stars, galaxies 
as well as in laser plasmas. The main instability condition of the several electrostatic and 
electromagnetic waves discussed here is the same as given in the form of inequality (24). 
Therefore, in our opinion, electrostatic perturbations are strongly coupled with magnetic 
field fluctuations in inhomogeneous un-magnetized plasmas. 
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